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Abstract: In this paper we consider piecewise affine differential equations mod- 
eling gene networks. We work with arbitrary decay rates, and under a local 
hypothesis expressed as an alignment condition of successive focal points. The 
interaction graph of the system may be rather complex (multiple intricate loops 
of any sign, multiple thresholds...). Our main result is an alternative theorem 
showing that, if a sequence of region is periodically visited by trajectories, then 
under our hypotheses, there exists either a unique stable periodic solution, or 
the origin attracts all trajectories in this sequence of regions. This result ex- 
tends greatly our previous work on a single negative feedback loop. We give 
several examples and simulations illustrating different cases. 
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Cycles limites dans les modeles afRnes par 
morceaux de reseaux genetiques incluant 
plusieurs boucles d'interaction 

Resume : Ce rapport porte sur des equations differentielles affines par mor- 
ceaux utilisees comme modeles de reseaux genetiques. Les taux de degradation 
sont arbitraires, et nous faisons unc hypothese locale d'alignement de points 
focaux successifs. Le graphe d'interaction du systeme pent etre assez complexe 
(multiples boucles d'interaction de tous signes, seuils multiples...). Le resultat 
principal est un tlicoreme d'alternative montrant que si une suite de regions 
est periodiquement visitee par les trajectoires du systeme, alors sous nos hy- 
potheses, soit il existe une unique solution periodique stable, soit I'origine est 
attractante pour toutes les trajectoires dans la suite dc regions. Cc resultat 
etend nos travaux precedents sur une boucle de retroaction negative. Plusieurs 
exemples et simulations sont presentes pour illustrer la methode. 

Mots-cles : systemes dynamiques lineaires par morceaux ; trajectoires pcr- 
iodiques ; applications monotones, concaves ; graphe d'interaction ; modeles de 
reseaux genetiques 
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1 Preliminaries 

It has been discovered in the 1960's that some proteins can regulate (i.e. ac- 
tivate or inhibit) the expression of genes in a Uving organism. Since proteins 
are the product of gene expression, feedback appears to be at the core of this 
process. Moreover, the unprecedented developments of cell biology in the last 
decades has led to a view where gene regulation involves huge numbers of ele- 
ments (genes, mRNAS, proteins ...), interacting in a nonlinear way. It is thus 
clear today that mathematical models and tools are required to analyse these 
complex systems. Several classes of models have been proposed to describe gene 
regulation. Although stochastic models undoubtedly have a major role to play, 
we will focus in this note on deterministic models. They can be divided in two 
classes : models based on differential equations [H], and discrete models based 
on a representation by a finite number of states [131 HHI ■ Both classes have their 
complementary advantages, the most obvious being relative to the accuracy vs. 
tractability dilemma. First proposed by L. Glass [4j, piecewise linear (in fact 
piecewise afSne) models appear to be an efficient intermediate between the two 
previous classes. Thus, the sometimes called 'Glass systems' have been both 
used to model real gene networks [3 [T7], and studied mathematically. The 
present paper is more related to this second type of work. 

The analysis of periodic solutions in piecewise affine gene network models is 
well characterised in the special case where all decay rates are supposed equal 
[SI [71 [51 [T5] . With non uniform decay rates however, very few results are known, 
and the techniques developed for homogeneous decay rates cannot be gener- 
alised. This lack of results is unfortunate, since real gene networks are known 
to involve very distinct degradation rates, both among mRNAs and proteins. 
Moreover, the question is also relevant from a more mathematical viewpoint, 
since systems with uniform decay rates display appreciably simpler dynamics 
than systems with distinct decays. Actually, in the uniform decay rate setting, 
trajectories are locally straight lines, whereas distinct decays lead to pieces of 
exponential curves. As an illustration of this, a 3-dimensional example of chaotic 
behaviour has been provided 114], whereas chaos in a uniform decay rates set- 
ting has been proved to require at least 4 dimensions [15j . 

In a previous paper [51 [TT] , we have successfully used tools from the theory of 
monotone systems and operators to tackle this problem. This has lead us to 
prove the existence of a unique stable limit cycle in an invariant region of state 
space that is periodically crossed by flow lines. This result holds only for a 
particular class of systems, namely those having a single negative loop as an 
interaction graph. 

However, the essential property of these systems that was mathematically used 
in the proof was a geometric condition. It was expressed as an alignment con- 
dition satisfied by pairs of successive focal points (see below for a precise def- 
inition), which was shown to be always true in the case when the interaction 
graph is a negative loop involving all variables of the system. In this paper, we 
generalise our previous results to a setting where the main assumption is only 
this alignment of successive focal points, regardless of the interaction structure. 
It is possible to show that this alignment assumption necessarily holds when 
each node in the interaction graph has outgoing degree 1. But the notion of 
interaction graph is essentially local, and the previous condition only has to 
hold at each point in state space (or in a region of interest in state space). This 
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includes of course feedback loops, but also many other interaction structures, 
including systems whose global interaction graph has a complex form, as some 
examples will illustrate. 

Since all these results require some technicalities, which have been partly ad- 
dressed in [Sm], the proofs provided in this paper are not entirely self-contained, 
but the missing parts appear in the cited references. 

The following section introduces the class of piecewise affine models of gene 
networks. Then in a next section, a general fixed point theorem for monotone 
systems is recalled, and applied to a class of piecewise affine systems. It is 
illustrated on several examples in a final section. 



2 Piecewise affine models 
2.1 Definitions and notations 

In this section we recall basic facts about piecewise affine models [H H [HI [2] . 
The general form of these models can be written as: 

^ = k{x) - Tx (1) 

The variables {xi...Xn) represent concentrations in proteins or mRNA pro- 
duced from n interacting genes. By abuse of language, the n elements of the 
system will be called genes in the sequel. Since gene transcriptional regulation 
is often considered to follow a steep sigmoid law, it has been suggested that 
idealised, discontinuous switches may be used instead to model these complex 
systems [4 . Accordingly, let us denote s+(- , 6*) : R {0, 1} the increasing step 
function: 

s+{x,e) = if2;<6', 
s+(a;,6i) = 1 if2;>6', 

which represent an effect of activation. Also, s~(a;, 6*) = 1 — s+(x,0), is its 
decreasing version, and represents inhibition. Unless further precision are given, 
we leave this function undefined at its threshold value 9. 

The map k : R" R!f: depends on x only via step functions of the form 
s^(a;i, 9i), and is thus piecewise constant. F G R"^" is a diagonal matrix whose 
diagonal entries Ta — 7^ , are degradation rates of variables in the system. 
As a concentration, each variable Xi is nonnegative and bounded. When Xi 
reaches a threshold value, this has an instantaneous effect on the system: for 
some j, the value of kj changes. In other words, a step function s^{xi,9i) 
appears in the expression of the function Hj. For each i G {1 • • • n} let us define 
a finite set of threshold values: 

e.-{0^...,C}, (2) 

where the thresholds are ordered: 9° ^ < 9l < ■ ■ ■ < 9^'^^ < 9f. The 
extreme values and 0^' are not thresholds but bounds on XiS value, but it will 
convenient to denote them as elements of Qi. 

Now, each axis of the state space can be partitioned into open segments between 
thresholds. Since the extreme values will not be crossed by the fiow (see later), 
the first and last segments include one of their endpoints : 

e [[91 9l), {9l}, {9l 9f), {9f},---A0r'}, (0T-\ 9f]] (3) 
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Each product D — YYi=i '^i defines a rectangular domain, whose dimension is 
the number of Di that are not singletons. When dimD = n, one usually says 
that it is a regulatory domain, or regular domain, and those domains with lower 
dimension are called switching domains [2] • We use the notation S! to represent 
the set of all domains of the form D above, and 2!r (resp. &s) for the set of all 
regulatory (resp. switching) domains. Since S^r is composed of finitely many 
domains, it will be convenient to identify it with 

n 

A = \{{Q---q,-l}. (4) 

1=1 

In the rest of the paper, this will lead us to use formulations like 'in a domain a'. 
This identification can be described by a map d : ^ yi, d {WM^'^ ^ C')) = 
(ai - 1 . . .a„ - 1). 

The dynamics on regular domains, or regular dynamics, is easy to describe, see 
next section. On sets of S!s on the other hand, the flow is in general not uniquely 
defined. This can be circumvented by using set-valued solutions and the theory 
of Filippov funs]. However, in this paper we will not need this theory, thanks 
to some mild assumptions explained in the next section. 

2.2 Regular dynamics 

On any regular domain D, the production rate k is constant, and thus equa- 
tion is affine. Its solution is explicitly known, for each coordinate i : 

^,{x, t) = x,{t) = ^ - e-^'* (x, - , (5) 

and is valid for all t G R+ such that x{t) e D. One sees above that the point 

V 71 In J 

It is an attractive equilibrium of the flow ([5]). Hence, if it lies inside D, it is an 
asymptotically stable steady state of system ([T]). Otherwise, the flow will reach 
the boundary dD in finite time, unless (/'(D) lies exactly on the boundary of T). 
However this situation is clearly not generic, and thus one will always suppose 
in the following that: 

When the flow reaches dT), the value of k (and thus, of changes, and the flow 
changes its direction. The point 4'{T>) is often called focal point of the domain 
T). Note that if d(©) — a, we will often denote it 0(a), or 0". 
It follows that the continuous trajectories are entirely characterised by their 
successive intersections with the boundaries of regular domains, and that this 
sequence depends essentially on the position of focal points. However, the def- 
inition of trajectories on the boundary of regular domains requires further ex- 
planations. Let us describe the case of singular domains of dimension exactly 
n — 1. Let be a n — 1 dimensional domain, intersecting the boundaries of 
two regular domains T> and D'. If in these two domains flow lines both point 
towards, or away from W (i.e. the flow coordinate in the direction normal to 
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W has different signs in D and CD') the latter is called respectively black wall or 
white wall. In both cases, the Filippov theory |12j . or some other technique, is 
required. Otherwise, i.e. when flow lines both cross W in the same direction, 
one usually call it a transparent wall, and trajectories on W can be defined by 
continuity, from D and D'. 

In the following, one will only deal with transparent walls. A simple criterion 
to ensure that all walls are transparent is the absence of auto-regulation, in the 
sense that no production term Ki depends on Xi. Then, the only regions that 
are excluded are the singular domains of co-dimension 2 or more, which form a 
rare set in state space (and even in f^^). 

Given the flow in a box D (of image a in yi), it is easy to compute the time 
and position at which it intersects the boundary of D, if ever. The position of 
the focal point with respect to thresholds determines entirely which walls can be 
reached: {x \ Xi — 6*"'"^} (resp. {x\xi = 0°*}) can be crossed if and only if (j^i < 
e^--^ (resp. > e^'). Then, let us denote /+ t(a) = {i e {1 • • • n}| 0, > 6*"'}, 
and 4-4(0) = {i G {1 . . < ^f'-^}. Similarly, /„„t(a) = /+,(a) U I'M 

is the set of escaping directions of D. Also, the following pairs of functions will 
be convenient notations: 0f:A^ 9^, 6'j^(a) = O^'^^ and 0^{a) — 0"'. 
When it is unambiguous, we will omit the dependence on a in the sequel. 
Now, in each direction i e lout the time at which x{t) encounters the corre- 
sponding hyperplane, for x € 'Da, can easily be shown to be: 

..(.). ^mfminl^^V^^^^^^ 

Then, t{x) — minig/^^^ ^i^x), is the time at which the boundary is crossed by the 
trajectory originated at x. Then, the escaping point of D from initial condition 
X takes the form ip{x, t{x)). Since this will be repeated along trajectories, x will 
generally lie on the boundary of the current box, except for the initial condition, 
which may however be supposed to lie on a wall without loss of generality. In 
this way, one defines a transition map T"" : dT) — > dT): 

T^x = tp{x,T{x)) , , 

= (j) + a{x){x-(j)). ^' 

where a{x) = exp(— T(a;)r). The latter depends on a, as seen from (j6|). 
The map above is defined locally, in a domain a. However, under our assump- 
tion that all considered walls are transparent, any wall can be unequivocally 
considered as escaping in one of the two regular domains in bounds, and incom- 
ing in the other. Hence, on any point of the interior of a transparent wall, there 
is no ambiguity on which a to chose in expression ([T]). In other words, there is 
a well defined global transition map on the union of transparent walls. Let us 
denote this map T. 

Now, the initial system ([T]) may been reduced to a discrete time dynamical sys- 
tem (DomT,T), where Dom^ is the subset of n — 1 domains of where all 
iterates of T are defined. From the previous discussions, it appears that Dom^ 
is the union of all (transparent) walls, minus the union of all finite-time preim- 
ages (i.e. finite number of backward iterates of T) of n — 2 dimensional (or less) 
singular domain. The topology of this domain is not trivial in general, and is 
described with more detail in 'W . 
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Now, we will focus on situations where there is a wall W , and a sequence . . .a^ 
of regular domains such that {T"-" o T"'"' • ■ • o T''''){W) V\W ^Qi. Then, fixed 
points of such an iterate are equivalent to periodic trajectories in the original 
system ((!]). 



3 An alternative theorem for convergence 

Let us state some general notations. We denote x < y and a; ^ y if these in- 
equalities hold for each coordinate (resp. entry) of vectors (resp. matrices) x 
and y. We call this order the partial order. Then, we denote x < y if x ^ y and 
X ^ y. For x ^ y, [x,y] — {z \ x ^ z ^ y}, and (x^y) ~ {z\x < z < y}. For any 
set A, by int {A) we denote the interior of A, and by A its closure. 



3.1 Theorems for monotone systems 

Let us state first the main theorem we want to apply. It is a fixed point theorem 
for monotone and concave operators, with respect to the partial order. Many 
variants of this theorem have been proposed since early works, more than 50 
years ago. The form we use is due to Smith jT^]. In words, it states that a 
monotone and concave map on a compact domain of may have either the 
origin as a unique fixed point, or a unique positive (for the partial order) and 
attracting fixed point. The second case happens when either the origin is an 
unstable fixed point, or when it is not a fixed point at all. Let us now state this 
theorem: 

Theorem 1. Let p G R", p> 0, andT : [0,p] — > [0,p] continuous, in (0,p). 
Suppose DT{0) — lim DT{x) exists. Assume: 

a:>0 

(M) DT{x) >0 if x>0, X <p. 
(C) DT{y) < DT{x) if < x < y < p. 
Assume also Tp < p. 

//TO = 0, let X = p{DT{0)), the spectral radius of DT{0). Then, 
A ^ 1 => \/x G [0,p], T^x — > when n oo. 

A > 1 => There exists a unique nonzero fixed point q ~ Tq. Moreover, 

(0,p) 

and for every x G [0,p\ \ {0}, T^^x q as n oo. 

Lf TO > 0, then T has a unique fixed point q G [0,p]. Moreover, q G (0,p) and 
T"a; ^ q as n ^ oo for every x G [0,p]. 

One may remark now that in the case when T has second-order derivatives, 
the concavity condition (C) admits a simple sufficient condition. Let the map 
Ti : [0,p] — » R+ denote the ith coordinate function of T : [0,p] — >• [0,p]. 
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Proposition 1. Suppose that for all z, j, /c G {1 • • ■ n}, and for all < x < p, 



dxkdx 



-(x) s; 0, 



and for all i,j there exists a k such that the inequality is strict. 
Then T satisfies condition [C] of theorem]^ 

Actually, it is clear that under this condition each term of DT is a de- 
creasing function of each coordinate x^. It is moreover strictly decreasing in at 
least one of these coordinates, and (C) thus follows. Observe by the way that 
the notion of concavity (w.r.t a partial order) we deal with here is weakened by 
the fact that it concerns only ordered pairs [x, y) of variables. 



3.2 Preliminary results for piecewise-afRne models 

Let C ~ {aP , ■ ■ ■ a^^^} denote a sequence of regular domains which is period- 
ically visited by the flow. Thus, we will consider — a^, and more generally 
the upperscript i in a* shall be understood modulo £, unless explicitly men- 
tioned. Also, we will denote walls in 6 as follows: — a^ D a*+i. Then for 
i € {0 • • •£} we use the notation 0* = 4){a'-), and 9^s} is the ith threshold, i.e. 

We define local and global transition maps as follows: 

rp(l) rp(2) rp(t) 

W" > > ■ ■ ■ > W*^ = W° (8) 

T=T<*'oT<'^"^>---oT<^> 

We denote a*+^ — a* = EiBg-, i.e. Si is the exiting direction of box a*, and 
Ei G { — 1,-1-1} indicates whether trajectories leave this box increasingly or de- 
creasingly in direction Bs- . 

Now, let us give a more detailed expression of a local transition map 
obtained after straightforward calculations: 

T«(x) = (<^} + (x,-</.})af(x))^^,. „ (9) 

where a^^\x) — a^^\xsi) 
is denoted T^'K 

Furthermore, we make the 

Assumption 1. For each z G {0 • • • ^ — 1}, for all J G {1 • • • ri} \ {si^i}, 0^^^ = 
4>'j . Pairs of successive focal points satisfying this condition will be said to be 
ahgned. 

^Remark that walls have be defined as closed sets. The transition maps are originally only 
defined on the interior of these regions. However, it is not difliicult to show that they can be 
extended at the boundary of these set in the present context. 




. The jth coordinate map of T*^*) 
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One may observe that there is a first dichotomy in theorem[T]: depending on 
whether is a fixed point or not, one has to check a condition on the spectral 
radius of DT(0) or not. Here, will be some corner of the wall we chose 
as a Poincare section, hence a point with all its coordinates equal to some 
threshold. As might be seen by some hand drawing, it appears that the case 
TO > happens when several parallel threshold hyperplanes are crossed in some 
direction, while is a fixed point if there is a single threshold in every direction 
(which is translated to 0). Let us state this fact more precisely. 

Proposition 2. Suppose that there are two distinct crossed thresholds in at least 
one direction. Assuming without loss of generality that one of these defines the 
wall i.e. sq is a direction with two crossed thresholds, one may write this 
condition as follows 

3i e {I- --n}, 3k £{!■■■£}, i = so = Sk and Of^^Of^. 

Then, rC^) o T('=-i) • ■ • o T^^) C int {W^) , and consequently T {W^) C 

int (W") . 

To ease the reading, we have postponed the proof of this proposition to An- 
nex [X] 



Remark 1. The previous proposition holds for any disposition of the focal 
points, but requires the assumption that each box in the sequence admits a unique 
escaping direction. 

Another important result that we need is the existence of a region on which 
T is monotone. If so, concavity will follow. Under the alignment condition 
on focal points, there will be a unique region on which this holds after a finite 
number of iterates. We shall prove this fact. 

First of all, let us compute the first order derivatives of a transition map T^*) 
at a point x: 

fa^*'' (x) if k = j 

otherwise. 
Where j G {1 • • - n} \ {s^-i} and fc G {1 • • - n} \ {sj. 

One deduces that diagonal terms of the Jacobian are positive, on the column Si 
one has 

( dT^^^ \ 

sign ( -Q^i^) 1 -sign ((/)^. - Xk) sign {(pl^ - Xs,) , fc S {1 • • • n} \ {s,}, 
and the Jacobian is zero elsewhere. 

This indicates a possible usefulness of partitioning each wall into zones of 
the form 



Z'(ct) = [x ^W^\ sign (0j - Xj) = cTj Vj} where a G {-1, +1} 
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And actually, we will show now that for a periodic sequence C = {aP ■ ■ ■ a^~^} 
with the two properties mentioned at the beginning (aligned focal points and 
one exit direction for each a*), and the fact that 6 is 'full dimensional', there is 
a single zone of interest on each wall. 

Proposition 3. Suppose that along C all variables switch at least once : {si\i £ 
{1 ■ ■ ■ i}} ~ {1 ■ ■ ■ n} . Then under the hypotheses above, for each i G {0 ■ ■ ■ £—1}, 
there exists a unique ct* £ such that 

• Va; G W\ T(*) o r(*-i) o • • •r(i) o • • • o T'^'+^\x) G Z'{a') : all orbits 
eventually enter this zone. 

• T^*) (Z*(cr*)) C Z*+^(cr*+^) ; no orbit escapes these zones. 

Proof. Let i G {!■ ■ ■£}. We are in fact going to define explicitly the sign vector 

First recall the alignment condition on focal points: 

yjj^si+u 0} = <^}+^ 

Recall also that Si defines whether is crossed increasingly or decreasingly in 
its normal direction, Sj. Then, one can then see from the definitions that: 

Si = sign {cPi^ - e^) = sign (<^',^ - a:,J, Mx G W'-\ 

Now, let x° G be an arbitrary point, and then define inductively x*"*"^ = 
T^^^x^, so that {x^}i is a trajectory in G. 

Let us introduce yet another notation: = sign(0J- — Xj). This quantity 

depends on .x°, at least at first sight. 

From the alignment condition and the expression of r(^+i) one gets: 

- = -"r''(^i+i) • -^))^ 

And since ot'^^^^ is nonncgative, we have Vj ^ Sj+i, A}"^^ = 

Furthermore, the expression of £j given above shows that in direction Sj, the 

precise value of coordinate x\. has no influence. In terms of our new notation: 

A|, = Si independently of x^, since by construction x\. = 6^s} ■ 
Then, one should remark that 

A^i = sign(</,^i-(</,^i+al:+i)(4^J.«-<^^i))) 

= sign - a;^) = sign (^^f - ^(f ) 

since a;' G C {a; | Xsi = O^s}}- Here the important point is that A^+^ does not 
depend on a;°, and an easy induction shows that this holds for all subsequent 
A™, when m > i + 1. 

So if all variables switch at least once, as we have supposed, A^ is independent 
of the initial condition, i.e. sign (T(a:°) — 0^) is fixed. Starting from other walls 
than does not change the argument. 

As for invariance, it is a consequence of invariance of the zones ^'(ct) under 
the action of the flow in each box, which can be retrieved from the explicit 
expression of this flow. □ 
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3.3 Main result 

We can now recapitulate the results of the previous section, and use them to 
apply theorem [T] This will be formulated as a single theorem, but before that 
we make some remarks. 

If the condition of proposition [3] is verified (all variables switch) , but not the 
condition of proposition [21 then it is not difficult to see that ^ single 

point, which is furthermore a fixed point of T. In the statement and proof of 
the theorem below, this point is denoted 0. 

When the condition of |)roposition[5]holds, some directions involve several dis- 
tinct thresholds, and flj^o ~ ^- -'^'^ ^^^^ case, will denote the corner point 
of W'^ (i.e. the boundary point with all its coordinates being threshold values), 
which also belongs to Z^{a^). 

Finally, it is clear that the map T is differentiable inside and that its differ- 
ential DT can always be extended by continuity to the point we conventionally 
denote 0. 

In words, the theorem below states that given a cycle of regular domains where 
successive pairs of focal points are aligned and all variables switch, there exists 
either a unique stable and attracting periodic orbit, or is the only attractor. 
The alternative depends on the stability of 0, and if furthermore two parallel 
thresholds are crossed, is not fixed and there is a unique stable periodic orbit. 

Theorem 2. Let C = {a", • ■ • a^^^} denote a sequence of regular domains 
which is periodically visited by the flow, and such that each domain a* has a 
unique exiting direction Si. Suppose that the focal points of Q satisfy Assumption 
[11 i.e. they are aligned. Suppose also that all variables are switching at least 
once. 

Consider the first return map T : W'^ defined in 13). Let X = p{DT{0)), 

the spectral radius of DT{0). Then, the following alternative holds: 

i) if X ^ 1, then Va; G T"^ when n oo. 

ii) if X > 1 then there exists a unique nonzero fixed point q ~ Tq. Moreover, 
q G int (Z'^(cr'^)) and for every x G \ {0}, T"x — > q as n ^ oo. 

If moreover the condition of proposition\^is satisfied, then the conclusion of 
ii) holds. 

Proof. To prove this statement, we verify that the hypotheses of theorem [T] are 
satisfied. First, from proposition [3l any x S enters 7l'{a^) under the action 
of T. Thus, we can consider the restriction of T to this zone from now on. We 
now give the main arguments for a proof that this restriction is monotone and 
concave, in the sense of conditions (M) and (C). 

From equation [ini it follows that the Jacobian of each local transition map has 
fixed sign in the zone Z'((T*), with nonzero terms only on the diagonal and on 
the column Si. From the chain rule, it follows that the Jacobian of T, denoted 
JT, is the product of matrices of this form. Since Si takes all values in {1 • • • n} 
by assumption, JT has no zero terms, and from proposition [3] one can easily 
deduce that these terms have a fixed sign on Z°(ct°). Then, by a simple change 
of coordinate system, it is always possible to ensure that these terms are in fact 
positive (see OUT] for an explicit example of this change of coordinates). 
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Now, the condition (C) can be proved to hold thanks to proposition [T] and the 
exphcit form of second order derivatives of local transition maps: 



(-) 



Ik Ofc (a:) 

7s, 4>s, - 



if k ^ j, m = Si or j 



7s, V 7si/ ((/-i, -XsiY 
otherwise. 

(11) 

where it appears that these quantities are of fixed sign in the zones Z*(cr*). Then 
we can apply proposition [U up to the same coordinate change as for the proof 
of (M). Because transition maps are monotone (up to coordinate change) it is 
possible to prove that the full return map T is also concave, just as the local 
maps. Here again, we refer to E] for a detailed justification. 

The statements i) and ii) are now obtained by strictly applying theorem 
[TJ the only missing hypothesis being Tp < p. The latter can be proved using 
similar arguments as those in the beginning of the proof of proposition [21 Here 
again, we refer to [TT] for a detailed proof. 

Finally, the last part of the theorem is a direct consequence of proposition [H 
Actually in this case, whatever corner point has been chosen as origin, the 
proposition shows that it is mapped in the interior of 2P{a^), which is identical 
to the condition TO > of theorem [1] whence the identical conclusion. □ 

This new result bears some resemblance with the previous work on so-called 
cyclic attractors [51 [B] . The important improvement is that we do not make any 
assumption on decay rates, which were uniform in this previous result. On the 
other hand, we have to suppose the alignment condition on focal points for our 
result to apply, whereas the cited references do not make any assumption on 
the precise position of focal points. 



4 Examples 

In this section we study three examples. They are not directly inspired by real 
biological systems, but serve the purpose of illustrating potential applications 
of theorem [21 



4.1 Two intricate negative loops 

Our previous result [51 [TT] can be seen as a particular case of theorem [51 for 
a system consisting of a single negative feedback loop involving all variables, 
there exists a cycle C satisfying the conditions of this theorem. Moreover in this 
case, we showed that i) holds with two variables (by showing A = 1), and ii) 
holds with three variables or more (by showing A > 1). 
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In this section we will consider systems with three variables consisting of 
two negative loops intricate in the following manner: 



First, the graph above corresponds in fact to various different systems. In 
particular, since gene 3 acts on both other genes, this may happen in general 
at two distinct threshold values, whose order needs to be specified. Moreover, 
the production rate of gene 3 is a function of X2 of the form s"'"(a;2), and the 
value Ks/js relative to these two thresholds can be chosen in two qualitatively 
distinct ways: between the two thresholds, or higher than the greatest one 
(excluding the case whore it is lowcir than the min of the two thresholds). Also, 
since gene 2 is regulated by the two other genes, we must also chose whether 
this happens as a sum or a product of step functions (a rapid inspection shows 
that in purely boolean terms, these are the only two cases in accordance with 
the graph above). In summary, there are thus eight cases to distinguish (two 
choices of thresholds orders, two choices of production rate and two choices 
of input function on X2)- 

In all cases, the set of regular domains is represented by ^ = {0, 1}^ x {0, 1, 2}. 
Any of the four choices fixes the position of the 12 focal points, and thus a 
discrete transition structure which may contain a cyclic sequence of the type 
studied in this paper. According to our inspection, in two cases this leads to 
a structure with no cycle, in one case to a structure with two cycles involving 
2 switching variables only, in four cases to a structure with one or more cycles 
presenting escaping edges (or walls) and finally in a single case to a structure 
with a cycle involving all three variables, and no escaping edge. In the following 
we analyse this last case only: it corresponds to 6l (resp. 9^) being the threshold 
of the activation 3 — > 1 (resp. the inhibition 3 H 2), and the regulation of X2 
being a product of step functions, with constraints on Ks as below: 

' Xi{t) = KiS+{x3,0l)--fiXi 

X2{t) = K2S-{xi)s-{X3,9l) -^2X2 

< X3{t) = K3S+{x2) -'jaxa 

with constraints Ki > Oi^i, i = 1,2, ^373 > K3 > 63 73. 



Up to a division by the degradation rates 7^ the focal points are given be 
the following table: 



e : 000 


001 


002 


010 


oil 


012 


100 


101 


102 


110 


111 


112 







Ki 







Ki 





Ki 


Ki 





Ki 


Ki 


K2 


K2 





K2 


K2 
































K3 


Ks 


K2 











K3 


K3 


K3 
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Figure 1: The limit cycle resulting (represented with a thicker line than the 
transient regime) from the example in section HTT] Parameter values are : 6i = 
02 = 01 = 1 and el = 2. (71,72,73) = (1,3,6). = 1.3, K2 = 5.1, = 10.8. 
Initial condition: x{Q) ~ (1.8,1.9,2.8) is chosen outside the region where the 
cycle lies, to indicate that the basin of attraction of this cycle is the whole space. 

Then, the transition structure mentioned above can be depicted as follows: 



110 ^111 ^ 112 




000 ^ 001 ^ 002 



where nodes represent regular domains, and arrows represent the transitions 
imposed by the flow (or the position of focal points). We have underlined the 
cycle C using bold arrows. It is easily seen that any of the 6 regular domains 
not belonging to C contain only initial conditions which enter 6 in finite time. 
In fact, we can remark that the 4 domains with 03 = 2, i.e. the half-space 
X3 > 6*1 is repelling for the flow. Hence, we can restrict the study of this system 
to the remaining half-space, which can be done be fixing s~(a;3,0|) = 1 in the 
equations. Then, we immediately see that the obtained system is a negative 
feedback loop with three variables. It is easy to check that the focal points in 
the cycle satisfy assumption [1] for instance from the table above. Thus, we 
can apply theorem [21 and moreover, we also know from [HUT] that the origin 
is unstable (i.e. A > 1), hence the conclusion ii) holds. In brief: the cycle 6 
contains a unique stable periodic orbit, which attracts all initial conditions in 
the 12 regular domains, as illustrated on figure [TJ 

4.2 More complex interaction graph 

In this example, we chose a more complex interaction structure, involving multi- 
ple loops of different signs. Namely, we apply our result to the following system: 
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X2{t) = K2 [s~(xi) s~(x3) + s+(a;i) s+(x2)] - 72 X2 

X3{t) = K3 [S~{xi) S+{X3) + S+{xi) S+(X2)] - 73 X3 



with constraints 



K, > 6*^74, i = 2,3. 

Ki<eau Ki>diju K'{>9i-fi. 



whose interaction graph is depicted on the left: its arrows correspond to the in- 
teractions appearing in the production term, so that there is no exphcit self 
loop on the node 1. No distinction between activation and repression is made 
in this graph, to avoid multiple arrows. Actually, the action of 2 on 1 can be of 
both types depending on X3 and similarly for several other arrows. 
The step functions act with a single threshold 9i for each variable Xi, not shown 
in the equations above. The eight parameters Ki and 7^ may take any value 
satisfying the consistency constraints. 

Since there is a single threshold per variable, the regular domains of this system 
can be represented by = {0, 1}'^. Up to a division by the degradation rates 
7i the focal points are given by the following table: 



e : 000 


010 


110 


111 


Oil 


001 


101 


100 


Ki 


Ki 


K[ 








A-' 




Ki 


K2 


K2 


K2 


K2 




















K3 


K3 


K3 


K3 









In a more geometrical way, this cycle can be depicted as: 

011-< 111 

001 ->-101 



-I- 



t 



010 1->-110 

000 -< 100 

The transitions between regular domains in ^1 = {0, 1}^. The dotted lines rep- 
resents a white walls: since no trajectory starting outside these walls can reach 
them, we can ignore them without difRculty. 



In the table above, we have disposed the domains in the order followed by 
any trajectory under the given parameter constraints. In particular, there is 
a cycle 6 involving all the regular domains of this system. Finally, it is easily 
seen in this table that two consecutive focal points only differ in the switching 
direction: they are aligned in the sense require by theorem [21 Moreover, the 
three directions are switching. Hence, our theorem applies and this system ad- 
mits a unique attractor, which may be either the 'origin' {Oi, 92,03) or a stable 
limit cycle. We have not been able to prove that the origin is unstable, but all 
numerical simulations we have performed with various parameter values have 
led to a limit cycle, as illustrated in figured 



Despite its complicated interaction graph, this system is such that on any wall, 
each variable modifies the value of at most one focal point coordinate. This 
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Figure 2: An example of numerical simulation of the example from section 
14.21 A limit cycle is observed, and represented with a thick line in the figure. 
Parameter values are : 0^ = 1 for all i. (71,72,73) = (1,3,6). Ki = 0.4, 
K[ = 1.02, K'{ = 2, K2= 6, K3 ^ 12. Initial condition: 2;(0) = (1.1, 1.1, 1.1). 

property is the origin of the alignment of focal point, and might be used to 
study more examples with a complex global interaction graph, which simplifies 
locally. 

Our last comment on this example concerns assumption [T] Actually, this align- 
ment assumption could be thought to imply that focal points are vertices of a 
rectangular parallelepiped. However, a rapid inspection of the table above will 
show to the reader that this is not the case in our example, and that it is diffi- 
cult to draw global conclusions from the local condition of focal point alignment. 



4.3 Multiple threshold values 

We provide now a simple two dimensional example involving two distinct thresh- 
olds in one direction. Namely, x\ can switch at two thresholds Q\ < 0^, and X2 
at a single threshold denoted 62 (and omitted in step functions involving X2)- 
This system writes 




xiit) = Ki [s-{xi,9l) s-{x2) + s+{xi,el) s+{x2)] + K[ s+{x^,e\) (X2) - 71 XI 

i2{t) = K2 [s+{xi,e\) S+{X2) + S+{xi,el) S-{X2)] + K'^ - l2 X2 

with constraints Ki > 6*^71, K[ > 6*^71, K'2 < 6*272, K2 + > 02^2- 

The aim of this model is to illustrate the last statement of theorem [5J First, 
let us show the production rate values (or focal point coordinates multiplied by 
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Figure 3: On the left: the transitions between regular domains in ^1 = {0, 1, 2} x 
{0, 1}. The dotted line represents a white wall: since no trajectory starting 
outside this wall can reach it, we can ignore this wall without difficulty. On 
the right: an example of numerical simulation of the example from section [4.21 
A limit cycle is observed. Parameter values are : 9\ = 62 — l^ and 9\ = 2 
(dotted hues represent these thresholds). (71,72) — (1,5). Ki — 1.5, K[ — 2.7, 
K'2 = 0.5, K2 = 6.5. 



decay rates) in a table, as in the previous example: 



e : 00 


10 


20 


21 


11 


01 
















K'2 


K'2 


K2 + K^ 


K2 + K^ 


K2 + K^ 


K'2 



Here again, a cycle C involving all regular domains exists for any parameter set 
satisfying the specified constraints. Moreover, any pair of successive focal points 
only differ in the switching direction, i.e. assumption [T] is verified. Hence, we 
may apply theorem [H and since 9l and Of are both crossed in C, we conclude 
that there exists a unique stable periodic orbit attracting all initial conditions. 
This fact is illustrated on figure [31 



We would like to stress the fact that this example has been chosen for its geo- 
metrical simplicity, which gives an easy intuition of our result. However, more 
complex cases could also be treated, and in particular we recall that the parallel 
thresholds may not be crossed successively. Let us illustrate this on a discrete 
transition structure, without entering into the detail of underlying differential 
equations. For example, if we now that all pairs of successive focal points are 
aligned in the bold cycle below: 



110 Ill >-112 

100 ■<-T 101 102 

t A I I I 

010 ->-011^ 012 

^ I Y 

000 ^ 001 -< 002 



then, we can conclude that there exists a unique stable periodic solution of the 
corresponding piecewise-linear differential equations, because two thresholds are 
crossed in the left-right direction. The above example may seem a little artificial, 
but it serves as an illustration of a property that holds in higher dimensional 
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systems, whose projections in 3 dimension might be at least as complicated as 
the diagram above. 

5 Conclusion 

We have presented and proved in this paper a theorem about periodic solutions 
of piecewise linear models of gene regulatory networks. This theorem relates 
discrete transitions between regular domains of these systems and their actual 
solutions. It does so under hypotheses of a local nature (alignment of pairs of 
successive focal points), hence allowing applications to a large variety of ex- 
amples, some of which have been presented in the last section. The alignment 
condition is related to the interaction structure of the system only locally, al- 
lowing complicated global interaction graphs to be handled within the present 
framework. 

A possible follow-up of this work would concern the control of gene networks, a 
topic arisen recently as a tool for synthetic biology. Actually, we have shown in 
a previous paper [10] that if production and degradation rates can be modified 
by an experimentalist (a fact modeled using input variables), then the prob- 
lem of forcing a discrete transition structure in piecewise affine models could 
be expressed as a simple linear programming (LP) problem. We applied this 
to the control of steady states, but it could also be apphed to the control of 
stable periodic orbits. Actually, if it is possible to control the discrete transition 
structure so that it presents a cycle, and if moreover we are able to impose that 
successive focal points are aligned (which amounts to adding constraints to the 
LP problem), then the result of this paper shows that the controlled system 
will have only two possible attractors: the origin or a stable limit cycle. This 
could be a powerful tool for regulating a system between two functioning modes. 

Appendix 

A Proof of proposition [2] 

Proof. Let us deal first with the more intuitive case when the two thresholds 
o'fa , ^i^-* are crossed successively, i.e. there are two consecutive parallel walls. 
In other terms, one assumes first that fc = 1. Then, for any x € W^, Xs^ — 
Xso = O^S 1 and eq. ^ gives 

T^'\x) = <i>] + af {x, - ^]), J ^ sr. 

It is easily checked that either (fso < oi]^ < or oi^J > O^]^ > (f)]^, since 
trajectories leave a" to enter a^. It follows that the scalar ct^j^'' (j^s!,"^^ £ (0,1). 

Since it does not depend on x, we abbreviate it into a^^^. Now, remember that 
si is by assumption the only exit direction from a^. This implies that for all 
j Si, denoting 0^ the thresholds bounding and in direction j, one 
has: 9j'>(l)j>9~, or equivalently 
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Multiplying by a'^^\ this implies 

Since the left-hand side above is the image of [9~ , 9^] by the translated map 
X Tj^\x) — 0], one deduces the expected inclusion: T^^^ (W") C int (W^^)- 

Now suppose k > 1. 
Since the walls are closed, connected sets, and all maps are continuous, to 
show that some wall is mapped in the interior of another wall, it is sufficient 
to show that none of its point is mapped on the boundary of the target wall. 
To achieve this, let us first prove that the following equality holds for any i £ 

{i-.-n: 

i 

r« o r(*"i) o • • • o r(i) {w") n dw = f]w'. (12) 

The proof is by induction. To initialise this induction, let us consider how the 
boundary of the target wall of the first transition map T*^^) is intersected. That 
is, we describe T'-^' (W^) fl dW^. Once again it will be useful to consider the 
translated map: 

Tj'\x)~<t>]^afUxs,)ix,-^]), j^si. (13) 

Let us denote 9 J < 9^ the bounding thresholds of in direction j. Note that 
these are also the bounding thresholds of W^, for j ^ sq- We also denote 9^^ < 
6*+ the bounding thresholds of in direction si , so that = IljXso [^7 ' ^ 
{^io^} X n,>J(^7,0t] and W = U,<sM,Ol] X {9^^} x U,>s^()7 ,0+]. 

Now, an image point encounters the boundary of if and only if Tj^\x) = 9^ 
for some j ^ si. Considering (fT5|) . this is equivalent to: 

T^'\x) - = 9f - d^] = ix,,) (x, - (14) 

for one of the two values of 9^. Now, observe that: 

• since si is the only exiting direction, one 

• aj^^ i^si) £ (0,1] for Xs-^ G [9j_^,9f_^], it is a monotone function of Xsi, and 
takes the value 1 only for Xsi = oi]^ ■ 

The first observation above implies that 9J — < < 9^' — for j 7^ si, and 
for any a G (0, 1) and Xj £ [9J ,9^'] this in turn gives: 

9- -cj,]<a- {9- - (/.]) ^ a ■ {x, - <j>]) a ■ {9+ - 4>]) < 9+ - 4>] 

Hence, one sees that (fH|) may only be satisfied if o'p {xsi) = 1, which from 

the second observation above occurs exactly for Xs^ = 6'ii''. But this defines the 
hyperplane bearing the target wall . In other words, the preimage of dW^ 
by T(i) is 

n 

w'c^w^^ X{[9-,9+] n 1 = 9f), xs, = 9i]^} . (15) 
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Moreover, from a 

W I'g/W j ^ 1^ the restriction of T^^) to W^^nVF^ is the identity. 
Hence, the set above is exactly T^^) n dW^. 

Remark that the case k = 1 treated previously could have been deduced from 
this property, since in this case one has parallel walls, and thus fl = 0. 
Under the assumption that we deal with boxes with a single outgoing direction, 
the argument above generalises to any local transition map as: 

t(*+i) [w) n dw'+^ = w'r\ w'+^. (i6) 

Now another useful observation is that 

W' n W'+^ = dW n dW'+^ (17) 
as may be seen from the explicit descriptions of each of these sets. 

Now, suppose that the induction statement holds for some i. Then, we 
have 

o t(*) o • • • o r(i) {w°) n dw'+^ cw'n w'+'^ 

by (HH) and the inclusion T^') o ■ • • o T^^) {W°) C W\ Then, because T^'+i) 
acts as the identity on W n the only points of T^*) o ■ • • o T^^) {W°) 

whose image by t'*"*"^^ intersects dW^'^^ must also lie in fl But 
this intersection is a subset of dW^ by (|17p . Hence, these points lie in fact in 
o ^.(j-i) Q . . . o 2^(1) (-p^o^ PI .^yj^j^;!^ ^jjg induction hypothesis equals 

fXj^Q Using again the fact that t(*+^) acts as the identity on this set, the 
expected T^'+i) o T^^) o • • • o T^^) (W^") n dW'+^ = fl^to follows. 

Now, it is not difficult to see that (fTS|) can be generalised to the above 
intersection as follows: 

k n 

f]w'^ l[[ej,e+] n {x \ x,. = ^ = o • • • fc} 

i=0 j=l 

and the main assumption of this proposition implies that this set is empty. □ 
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